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CLUSTER AND VIRIAL EXPANSIONS FOR THE 
MULTI-SPECIES TONKS GAS 

S. JANSEN 


Abstract. We consider a mixture of non-overlapping rods of different lengths 
moving in M or Z. Our main result are necessary and sufficient convergence 
criteria for the expansion of the pressure in terms of the activities Zk and the 
densities pk- This provides an explicit example against which to test known 
cluster expansion criteria, and illustrates that for non-negative interactions, 
the virial expansion can converge in a domain much larger than the activity 
expansion. In addition, we give explicit formulas that generalize the well- 
known relation between non-overlapping rods and labelled rooted trees. We 
also prove that for certain choices of the activities, the system can undergo a 
condensation transition akin to that of the zero-range process. The key tool 
is a fixed point equation for the pressure. 

Keywords: Cluster and virial expansions, non-overlapping rods, combinatorics 
of labelled colored trees, close-packing transition. 


1. Introduction 

The present article deals with one-dimensional systems of non-overlapping rods 
on a continuous segment [0, L\ or a discrete interval {0,1,..., L — 1}. There are 
countably many types k of rods, coming each with a length If. > 0 and an activity z^- 
A rich literature deals with related models: our model is a multi-species variant of 
the well-known Tonks gas |Ton36] . We may also view it as a one-dimensional special 
case of hard spheres mixtures I.IKi 1 . A good control of discrete one-dimensional 
partition functions enters as a building block for two-dimensional models with ori¬ 
entational long range order [IVZ06[ IDU13j . The one-dimensional discrete system 
of rods also appears in stationary distributions for driven one-dimensional sys¬ 
tems [KLM + 02] , which in turn are closely related to the zero-range process where 
particles are piled up rather than aligned in a rod )EH05| . Phase transitions for 
one-dimensional cluster models have been studied in detail by Fisher and Felder- 
hof IFF701 lFis72l . 

Our principal motivation comes from the model’s solvability and the specific, 
though model-dependent, answers to questions on cluster expansions it allows. 
The first question concerns domains of convergence. There are many sufficient 
convergence criteria available, but it is an ongoing effort to improve them, see for 
example [FP071 IBFPlOj . This raises the question of how much room for improve¬ 
ment there actually is. We answer this question for the multi-species Tonks gas by 
determining the exact domain of convergence (see Eqs. 0 ctnd. (|2|) below). The 
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answer for general models can of course be quite different, but we hope that our 
results will serve as a helpful control group in future studies. 

The second question is how the domain of convergence of the activity expansion 
compares to the domain of analyticity. It is common wisdom that for repulsive 
interactions, the activity expansion ceases to converge before a phase transition 
occurs [R.uc69j . For single-species model this means that the radius of convergence 
is strictly smaller than the activity value at which a phase transition occurs, if it 
occurs at all. We prove that for the multi-species Tonks gas with rod lengths £k = k, 
the difference between convergence and analyticity domain is even more drastic: for 
the convergence of the activity expansion it is necessary that the activities go to 
zero exponentially fast Zk = 0(exp(—ak)) —> 0, while the pressure stays analytic 
all the way up to exponentially diverging activities Zk —> oo (Corollary 12.711 . 

The third question concerns the virial expansion, i.e., the expansion of the pres¬ 
sure in terms of the densities pk■ It has been suggested that the virial expansion 
can be more advantageous than the activity expansion |Bryll| , and indeed for some 
models this is known to be true [Joy88j [BMHj ■ We prove that the same holds for 
the multi-species Tonks gas, again in a quite drastic way: the virial expansion 
converges in all of the analyticity region lTheorems l2.5l and l2.121 Corollary 12.711 . in¬ 
cluding densities pk{zi, 22 , • ■ •) that correspond to exponentially diverging activities 
Zk —> 00 (when £k = k). 

In addition, we provide explicit formulas for the pressure-activity expansion that 
are interesting from a combinatorics point of view (Theorems 12.41 and 12.1111 . In the 
continuous case, we find that the activity expansion is (up to signs) the multivariate 
exponential generating function for labelled rooted colored trees. The vertex colors 
correspond to rod types and the trees have weights that depend on the lengths £k- 
This generalizes the well-known relationship between the exponential generating 
function T(z) = X)n>i z n n n ~ l /n\ of rooted labelled trees and the pressure for non¬ 
overlapping rods of length 1 (see IBI03a' and the references therein). It would be 
interesting to know whether the answers to corresponding combinatorial puzzles 
given in [BerOSl ITatlbj extend to the multi-species setting. 

Let us describe in more detail our results on the convergence domain of the 
activity expansion of the pressure in the infinite volume limit. In the continuous 
case, the expansion converges absolutely if and only if ('Theorem 12. 4l PI 

OO 

3a > 0 : ^ \zk\e aIk < a. (1) 

k =1 

In the discrete case we may choose ik — k, and the cluster expansion for the pressure 
converges if and only if (Theorem 12.111) 

OO 

3a > 0 : J2 N e ° fc < e“ - 1- (2) 

fc=1 

The principal novelty lies in the only if part: if the criterion fails, then the expansion 
is not absolutely convergent. Eqs. m and (0 should be compared with the following 
known sufficient criteria. In the continuous case, it is enough that for some a > 0 


^The criterion 0 refers to the pressure-activity expansion p(z). The convergence of the 
density-activity expansion pk(z) in general requires a strict inequality Yk \ z k\ ex P( a ^fc) < a • The 
same remark applies to the discrete system. 
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and all k e N jPTTMl 

OO r. OO 00 

Y / |zj|e a ^l([a;,a; + 4 ] D [0,4] 4 0]dx = ^(4 +4)^e“^ < o4- (3) 

In the discrete case, it is enough that for some a > 0 IGK71I lFP07ll 

OO 

Y |gfc|l({g, ■■ • ,X + fc - 1 } 5 0)e afc = y^fce afc |z fc | < e a - 1 . (4) 

x£Z,, k £N k —1 

The activity domains determined by the sufficient criteria ([3]) and are clearly 
smaller than the full convergence domains given by © and ©, however we shall 
take the point of view that the difference is small: in the discrete case both © 
and ® are of the form ||z|| a < exp(a) — 1 with weighted norms ||z|| Q that impose 
the same exponential decay of Zk and differ merely by the prefactor k. From 
this point of view the one-dimensional model of non-overlapping rods provides an 
example for which the classical convergence criteria are already nearly optimal 0 

A full list of results is given in the next section. In addition to exact formulas 
for the activity and density expansions as well as their domains of convergence, 
we prove that the pressure solves a fixed point equation. In the continuous case it 
reads 

P(z) = Y Zk exp (-4 p(z)^ (5) 

k 

and is satisfied by the pressure whenever the equation has a solution f Theorem 1 2.2 II . 
In the discrete case the equation is instead (Theorem 12.101) 

1 - exp(—p(z)) =Y Zk exp {-kp(z)). (6) 

k 

(remember £k = k). Situations where the fixed point equation has no solution are 
possible and correspond to a close-packing regime (Theorem 12. .'ll) . Section OH pro¬ 
vides two different explanations of the fixed point equation, a statistical mechanics 
explanation and a probabilistic explanation in terms of renewal equations. 

The fixed point equation is not only of interest in itself but also lies at the heart of 
all of our results and their proofs; it is no coincidence that Eqs. © and © resemble 
so closely the convergence criteria © and © as well as functional equations for 
combinatorial generating functions. This will become clear in the proofs, given in 
Sections Q] to [7] Section [5] proves a technical result (Theorem 12.81) on the inversion 
of the density-activity relation via the inverse function theorem. 

2. Model and results 

2.1. Continuous system. Let (4)fceN be a family of rod lengths 4 > 0 and 
z = (zfe)feg n positive activities Zk > 0. The rod lengths do not need to go to 
infinity; in fact the case 4 —> 0 might be of interest in view of the transition 
towards dust studied in some fragmentation models iHaa04] . Let I C Nq be the 
set of multi-indices (JVk)keN that have none or finitely many non-vanishing entries, 
and I* := fZT\{0} the set of multi-indices with at least one non-zero entry. For 

2 This point of view focuses on qualitative features of the domain of convergence for objects of 
unbounded size (L : —> 00 ). A different question is about quantitative estimates on the radius of 
convergence for hard-sphere systems of bounded size, see e.g. IFPS07I . 
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nelwe set n! := J} fc (nfc!) an d z ” := Tlfc z k k i with the convention 0 ° = 1 and 
0! = 1. The grand-canonical partition function for the volume [0,L] is 

Z L (z) := 1+ V' — — / dxn---dxijv 1 dx 2 i---dx 2 jv 2 --- 

m J[0,L]Sfe N k 

x l(yk,j : [x k j,x k j + 4] c [0, L]^ 

X ^ (?7l, z) I [x/jjf, Xfcj -f* 4] Pi H” 4i] = 0^ • (7) 

In the integral Xfcj is the left end point of a rod of length 4- The indicators ensure 
that rods do not overlap and lie entirely in [0,L]. The pressure is 

p(z) := lim y log S L (z) ( 8 ) 

L—too 1 j 

The existence of the limit in [0, oo) U { 00 } follows from general subadditivity argu¬ 
ments, a more concrete determination is given in Theorem 12.21 below. We shall see 
that the stability condition 

39 € R : Zk exp($4) < 00 (9) 

k 

ensures that the limit defining p(z) is finite. Note that we allow for 9 < 0; in 
particular, for 4 = k, the activities are allowed to diverge as exp(cfc). Let 9* 

9* := sup{0 G K. | Zk exp($4) < 00 } (10) 

k 

be the abscissa of convergence of the Dirichlet type series 

g{0) = ^2 z k exp(04). (11) 

k 

For notational convenience we suppress the ^-dependence in 9* and g(9). Depending 
on the values of z, the fixed point equation g[9) = —9 may or may not have a 
solution. As we shall see, the cases correspond to different physical behaviors and 
the following names are convenient. 

Definition 2.1. Let z = (zk)keN be non-negative activities that satisfy the stability 
condition ©. Then z belongs to the 

(a) fluid domain Pfl u id if 9* = 00 or 9* < 00 and g(9*) > —0*; 

(b) close-packing domain if 9* < 00 and g(6*) < — 9*. 

(c) transition domain if g(6*) = — 9*. 

In the fluid domain the fixed point equation g(9) = — 9 has a unique solution 
9 < 9* , in the close-packing domain it has no solution, and in the transitional 
domain the unique solution is 9 = 9*. 

Theorem 2.2. Let (zk)k^N be non-negative activities satisfying the stability con¬ 
dition © . Then 

(a) In the fluid domain the pressure p(z) is given by the unique solution p > —9* 
of z k exp(-p4) = p . 

(b) In the close-packing and transition domains the pressure is p(z) = — 9 *. 
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The theorem is proven in Section |4j an explanation of the fixed point equation 
is given in Section [3] Theorem 12.21 suggests the possibility of phase transitions, 
illustrated by the following example. 

Example. Take rod lengths Ik = k and parameter-dependent activities Zk{p) = 
exp (kp — Vk), /ifR. The activities are associated with a parameter-dependent 
Dirichlet series ^ fe ex P[(^ + p)k — Vk\- The abscissa of convergence is 8* = —p 
and at —8* the Dirichlet series takes the /^-independent value ^ fc exp(— Vk). In 
the fluid regime p < Y2k exp(— Vk) the pressure p(p) is the unique solution to 
p = Y2k exp (kp — kp — Vk), and an implicit function theorem shows that p(p) is 
analytic with 

dp = X; fc fcexp(% -p] - Vk) < Ylk kex P(-V^) 1 

d/i 1 + J2k kexp(k[p - p] - Vk) ~ 1 + J2k kexp(-Vk) 

In the close-packing regime p > Y2k exj P(~Vk), the pressure is p(p) = p, and we 
recognize a first-order phase transition at p = Y2k exp (—Vk). 

Remark. The previous example is easily extended. Let (ak)ke^ be non-negative 
weights such that h(z) = akZ k has positive radius of convergence R > 0. Set 
Zk{p) '■= a,k exp (kp). Then there is a phase transition in p if and only if R is finite 
and Y2k a kR k < oo. The transition takes place at exp (p) = R. If Y2k kakR k < oo, 
it is of first order. The situation is closely related to condensation in the zero- 
range process |EH05j . with the important difference, however, that we can take 
exp(/z) > R and go beyond the coexistence region, into the close-packed phase. 
The phenomenon is also similar to phase transitions studied in Fisher-Felderhof 
clusters |FF701 lFis72| . 

A very heuristic explanation of the phase transition is the following. Suppose 
that Ik = k and the abscissa of convergence 8* = — limsup^^^ ^ log Zk is finite. 
Then exp (-L8*) is the weight of the configuration where space is filled with one 
long rod. In order to see how much we loose by breaking space into smaller rods, it 
seems natural to rescale the activities as Zk —► Zk exp (k8*). If we neglect excluded 
volumes, the partition function for a system of many rods for the rescaled activities 
becomes exp(L YV°=i Zk exp(fc0*)) = exp (Lg(8*)). When g(8*) < —8*, it seems 
more advantageous to fill space with one long rod, and so we recover the case 
distinction from Theorem l2.2l though the use of rescaled activities in this argument 
is somewhat ad hoc. 

The next theorem shows that the fluid and close-packing domains do indeed 
correspond to different behaviors of the system. The relevant order parameter is 
the packing fraction, the fraction of volume covered by rods. The convergence 
Y2k^ k ^ r k/L —;► cr stands for convergence in probability in the grand-canonical en¬ 
semble, i.e., for all £ > 0, the grand-canonical probability that | J2k Nkik/L — a\ > e 
goes to 0 as L —> oo. 


Theorem 2.3. Let (-Jfc)fceN be non-negative activities satisfying the stability as¬ 
sumption ©• Then as L oo the packing fraction behaves as follows: 

(a) In the fluid regime 


Y yVfcjVfc u(z) < 1, <t(z) 

k 


Y2 k 4£fc exp[-4p(z)] 

1 + Efc ^ kZk exp[-4p(z)]' 
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(b) In the close-packing regime, the packing fraction converges to 1. 

The theorem is proven in Section [4] by a large deviations approach. 

Remark. In the transition regime g(9*) = —9 * there are two possible scenarios: if 
as 9 approaches the abscissa of convergence the derivative g'{6) = Y k Ik exp(#£fc) 
diverges, the packing fraction converges to 1. If the derivative stays bounded, let 
cr* := limg g'(9)/[\-\-g'{9)\. Then cr* £ (0, 1) and Lemmas l4.2l and l4.3l below only 
show that the grand-canonical probability of seeing a packing fraction smaller than 
a* — e goes to zero, for every e > 0. Intuitively, this corresponds to a coexistence 
region between a densely packed phase (cr = 1) and a fluid that saturates at cr = cr*. 

Now we give the domain of convergence of the cluster expansion, as discussed in 
the introduction, and an explicit formula for the expansion. The formula generalizes 
a well-known relationship between a tree generating function and the pressure for 
non-overlapping rods of length 1, see |BI03aj and the references therein. 

Theorem 2.4. Let (zk)ke n be non-negative activities satisfying the stability as¬ 
sumption ([9]) . The following holds: 

(a) IfYk z k exp(a4) < a for some a > 0, then the pressure defined by Eq. © 
is given by 

p(z)= E —( 12 ) 

z —' n\ z —' 

n.GX* k 

and the sum is absolutely convergent. 

(b) If Y k z k exp (a£k) > a f or a H a > 0, then 

n€l* k 

The formula (1121) is proven in Section [5j where we also provide a combinatorial 
interpretation in terms of tree generating functions. The convergence is addressed in 
Section [6] We should stress that (a) refers only to the convergence of the expansion 
of the pressure; for the convergence of the expansions of the densities p k {z), the 
inequality m in general has to be strict for some a > 0. 

Remark. For negative (or complex) activities with Yk l z fc| ex P( a ^fc) — a f° r some 
a > 0, the sum (fl2l) is absolutely convergent as well and we adopt Eq. 03 as a 
definition of the pressure. In Section[5]we will see that F(z) = —p(—z) satisfies the 
flipped fixed point equation F = Yk z k exp(^fcF); in fact the convergence condi¬ 
tion m is equivalent to the existence of a solution to the flipped equation. This is 
interesting because F is the pressure for a three-dimensional multi-type branched 
polymer |BI03b . 

Finally we examine the virial expansion, i.e., the expansion in terms of the 
densities 

p k (z) ■.= Zk-^-(z). 

OZk 

Theorem 2.5. Let z = (zfc)fcgN be non-negative activities in the fluid domain. 
Then the partial derivatives -§^{z), k £ N, at z exist and we have 

( ■, = J2 k Pk(z) „ = Pk(z) exp{£ k p{z)) 

VV) l-YklkPk{z)' k l-YjtjPi(z) 
with Y k ^kPk(z)) = <j{z) < 1 as in Theorem \2.d]f a). 


( 13 ) 
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Theorem 12.51 is proven in Section [7] by applying an implicit function theorem to 
the fixed point equation (0 • This approach also shows that the pressure is analytic 
in all of 2?fl u id- 

Corollary 2.6. Fix non-negative activities z° = (zjt)kcN £ £Wid- Then for every 
k £ N, the map Zk K > p(zi, • • •, zj*._ 1: Zk,z%,...) is analytic in some neighborhood 

of 4 

A similar statement holds for suitable parameter-dependent activities. In par¬ 
ticular, the parametrization Zk(t) = tZk allows us to connect vanishing activities 
to every z £ 2?fl u id in such a way that the pressure is an analytic function of the 
parameter t. 

Proof. For simplicity we consider k = 1, the other cases are similar. Let 9* be 
the abscissa of convergence of z% exp (4$) and f(z\,0) := 9 + z±exp(9li) + 
J2k >2 ~fc exp (9(h)- Then / is a holomorphic function of (zi, 9) in Cx {9 | Re# < 9*}. 
We have f(z®,—p(z 0 )) = 0 and 

%(*?> =1 + J2h4e~ thp{x0) = 1 _l (z o ) ± °- 

The holomorphic implicit function theorem [FG021 Chapter 7] (see also [Sok09j for 
a power series approach) guarantees the existence of a holomorphic function 9(z\) 
defined in some open complex neighborhood of 0 such that 0(zi) = — p(z° ) and 
f(9(zi), z\) = 0. When z i is real and close enough to z±, then z = (zi, zjj, z °,...) 
must be in TWid- Indeed, changing only one coefficient z\ does not change the 
abscissa of convergence 9 *, and the inequality z k exp(4#*) > —#* holds by 
continuity when z ° is replaced by z\ sufficiently close to 9*. It follows that p(z) 
solves the fixed point equation and therefore z\ H > p{z\, z$, z® ...) = —Q{z\) is 
analytic. □ 

It follows that the domain of convergence of the activity expansion is in general 
smaller than the domain of analyticity 2?fl u id- In contrast, the pressure-density ex¬ 
pansion converges absolutely when Yk 4#fc( z ) < 1> which is the case in all of T’fluid- 
As mentioned in the introduction, the difference between the domains is quite dras¬ 
tic: when Ik = k, the fluid domain and the activities for which the virial expansion 
converges can include diverging activities Zk —> oo, while the convergence of the 
activity expansion requires exponentially decreasing activities Zk = 0(exp(— ak)). 
Precisely, let 

X> May := {z £ R+ | 3a > 0 : ^ Zk exp(fca) < a} 

k 

X^vir := {z £ | z £ X>fl u id, Pk{z) exists for all k, and k Pk(z ) < !} 

k 

be the domains of convergence for the activity and the density expansions when 
4 = k. 

Corollary 2.7. Assume 4 = k for all k £ N. Then 

X^May ^ V vir — ^fluid- 
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Proof. Let pk = c/k 3 with c > 0 small enough so that J2k k P k < 1- Let p := 
(Efc Pk )/(1 -J2k k Pk) and z k ■= Pk exp (kp)/{l-J2j 3Pi)- Thenp = J2 k Zk exp (-kp) 
and it follows that p = p(z) and pk = Pk{z). Since ^2 k kpk < 1 , the virial expansion 
converges and z is in T> v i r . At the same time Zk diverges as k —> oo exponentially 
fast, so the activity expansion cannot converge and 2 ^ T> Ma , y . Thus V M ay jj) 2? v i r - 
The identity = Pfluid is part of Theorem 12.51 □ 

We conclude with a technical result concerning the use of inverse function the¬ 
orems which complements Section 2.2 in [JTTU14] . Often virial expansions are 
obtained by inverting the density-activity relation with the help of an inverse func¬ 
tion theorem. This works for finitely many species because dpk/dzj at 2 = 0 is 
the identity matrix and in particular invertible. For infinitely many species, it is 
still true that the matrix ( dpk/dzj ) of directional (Gateaux) derivatives at 2 = 0 
is the identity matrix. However the existence of directional derivatives is no longer 
enough to guarantee Frechet differentiability in suitable Banach spaces. In theory 
one could imagine that this is just a technicality requiring additional estimates. The 
next theorem shows that to the contrary, for the natural Banach spaces at hand, 
the usual inversion procedure cannot work because z i—» p(z) is not a bijection 
between neighborhoods of the origin. 

We restrict to £k = k. For a > 0, let ||a;|| a := \ x n\ exp(an) and E a := 

{(iCnjngN 6 R N | IMU < oo}. These Banach spaces are natural because of the 
convergence criterion || 2 || a < a from Theorem 12.41 They allow for negative and 
complex activities and densities. 

Theorem 2.8. Suppose that ik = k for all k £ N, and let E a ,Vb C C N be the 
spaces of complex activities and densities introduced above. There is no way to 
choose a > 0, b > 0 and neighborhoods U a C E a and V), C Eb of the origin so that 
p(-) is a bijection from U a onto Vj,. 

Theorem l2.8l is proven in Section [HI It explains why the virial expansion is much 
more delicate for infinitely many species than for finitely many species. A general 
result on multi-species virial expansions was nevertheless proven in [JTTU14] using 
Lagrange-Good inversion IGoofiOl IEM94I . 

2.2. Discrete system. On the lattice we may take without loss of generality £k = 
k, k € N. We use the multi-index notation from the previous section. The grand- 
canonical partition function for rods on a line {0,1,..., L — 1} with L £ N is 
defined as in Eq. 0, except that the integration over rod end points is replaced by 
a summation over end points in {0,..., L — 1}, and the indicators should ensure that 
rods are contained in {0,1,..., L — 1}. The definition of the pressure, the Dirichlet 
type series g{9), the stability condition ©, and the abscissa of convergence (fTOl) 
are unchanged. The fixed point equation changes, however, and we define 

f(0) = - log(l - g(0)) = - log(l - ^2 Z k exp(fc0)) • 

k 

By expanding the logarithm, we find that f(9) is a power series in exp (9) with 
positive coefficients and radius of convergence R < exp (9*). The equation f[6) = 
—0 has a solution 9 < log R if and only if g{9*) > 1 — exp (9*). 

Definition 2.9. Let (zk)k &n be non-negative activities satisfying the stability con¬ 
dition |9|). We say that the discrete system of rods is in the 
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(a) fluid regime if g(9*) > 1 — exp(0*); 

(b) close-packing regime if g(6*) < 1 — exp(0*); 

(c) transition regime if g(9*) = 1 — exp (9*). 

Theorem 2.10. Let (zk)ke n be non-negative activities satisfying the stability con¬ 
dition © • Then 

(a) in the fluid regime the pressure p(z) is the unique solution of f(—p) = p, 
and we have p{z) > — 9 *. 

(b) In the close-packing and transition regimes the pressure is p(z) = —9 *. 
The theorem is proven in Section [I] 

Theorem 2.11. Let (zk)ke n be non-negative activities satisfying the stability con¬ 
dition ©. The following holds: 

(a) If Yk z k exp(afc) < exp(a) — 1 for some a > 0, then the pressure defined by 
Eq. © is given by 



(14) 


and the sum is absolutely convergent. 

(b) If Yk Zk exp(afc) > exp(a) — 1 for all a > 0, then 



The explicit formula for the expansion coefficients is proven in Section [5j con¬ 
vergence is addressed in Section [G] 

Remark. There is an interesting activity region intermediate between the domain of 
convergence and the fluid domain. Write X and Y for discrete intervals {x,... ,x + 
k — 1}, {y,..., y + j — 1}. Suppose that 



Then the infinite volume Gibbs measure can be constructed as the unique stationary 
measure of a Markov birth and death process, the loss network jFFGOl] , The one¬ 
dimensional model confirms that the loss network representation works in a domain 
that is in general larger than the domain of convergence of the activity expansions, 
since the loss network does not require exponentially decreasing activities Zk = 
0(exp(— ak)). 

Theorem 2.12. Let z = (zk)keN be non-negative activities in the fluid domain. 
Then the partial derivatives -§§-(z), k £ N, at z exist and we have 



Pk(z) exp [(A; - 1 )p(z)\ 

1 ~Yk k Pk( z ) 


with Yk k Pk{z) < 1. 

The proof is analogous to the continuous case treated in Section [7] and there¬ 
fore omitted. A heuristic derivation of the pressure-density and density-activity 
relations is given in Section (3© 
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The pressure-density expansion converges again in all of 2?fl u id- In the special 
case when there are only monomers, i.e., Zk = 0 for k > 2, we recover the well-known 
equations p = log(l + z\), zi = pi /(1 - pi). 


3. A FIXED POINT EQUATION 


In this section we provide two different explanations of the fixed point equations 
for the pressure, a physical explanation in terms of van der Waals mixtures and a 
probabilistic explanation in terms of renewal processes. 


3.1. Van der Waals mixtures. We assume the reader is familiar with notions 
from statistical mechanics such as the constant pressure ensemble, and employ 
common approximations such as log A! ss N(\ogN — 1); for notational simplicity 
we follow physics conventions and pretend that the approximations are identities. 

Consider first the continuous setting. The partition function of the constant 
pressure ensemble at fixed number Nk of rods of size k can be computed as 

Q(N 1 ,N 2 ,...-p)= f ^ Je-pEAAjp-Fdr)"' 1 (15) 

where M = Nk is the total number of rods. The multinomial coefficient rep¬ 
resents the number of ways to assign rod types to the M rods labelled from left 
to right. The integral f n exp(—pr)dr comes from integrating over the spacing 
between two consecutive rods. Eq. (fT5l) yields the Gibbs free energy 

G(N U N 2 , ... ;p) = Nk (log + Ph) ■ (16) 

fc 

Alternatively, we can compute directly the Helmholtz free energy 

FW,iV 2 ,... V) = (log v _ £ - l). (17) 


Isee Lemma l4.ll) . Recall G 


F + pV with V = — dG/dp , p 


p = 


V-Ek^Nk' 


—dF/dV, which give 
(18) 


This is a multi-species variant of the van der Waals equation, and generalizes the 
well-known equation of state of the single-species Tonks gas |Ton36j . Using the 
relation log,?*, = dG/dNk = dF/dNk , we obtain 

Zk = ex P(P4)- (19) 

Eqs. m and (g9D explain the formulas in Theorem l2.5l The fixed point equation d5j) 
is obtained from Eq. (1191) by multiplying both sides with exp(— p(-k) and summing 
over k. 

For the discrete system, the integral in Eq. m has to be replaced by a geometric 
sum. The Gibbs energy (ITCl) becomes 


G(N i, N 2 ,... ;p) = ^ N k (log 


(1 - exp(— p))N k 

E,- Nj 



The equation of state is 


E k ^k 
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and the density-activity relation is 

(1 - exp(-p))N k 
z k = - ^-r - exp (jpk), 

compare Theorem 12.121 We multiply with exp(— pk), sum over k , and obtain the 
fixed point equation ([6]). 

3.2. Renewal theory. As done in [TVZOfi] . one-dimensional polymer partition 
functions can be treated by renewal theory [Fcl71 . Chapter XI]. The key idea is 
to reinterpret the line as a time axis and the starting point of a rod as an event 
(a light bulb breaks and has to be renewed), and the intervals between two events 
as waiting or interrarival times. The Gibbs measure is invariant with respect to a 
suitable rescaling of the activities IGK711 . If the activities can be rescaled in such 
a way that they define a probability measure on waiting times (the light bulb’s 
lifetime), then the canonical partition function is identified with the probability of 
a certain event and probabilistic techniques apply. The crucial point now is that 
this rescaling is possible if and only if the fixed point equation has a solution. 

For details, consider first the discrete case. For X C Z set <f>(Xr) = Zk if X = 
{x, x + 1,..., x + k — 1} is a discrete interval of cardinality k > 2, $(X) = 1 + zi 
if X has cardinality 1, and $(A r ) =0 otherwise. Then 

s L (z)= y *m---*(x D ) 

where the sum runs over all partitions of {0,1,..., L — 1}. An element Xj = {x} of 
cardinality 1 corresponds to either an unoccupied lattice site or a site occupied by 
a rod {x} of cardinality 1. Rescaling the activities as <f>£(AT) = £# x <fr(X) for some 
£ > 0 multiplies the grand-canonical partition function Eb^z) by £ L and leaves the 
associated probability measure unchanged |GK71 i. We wish to choose £ so that 

OO 

(1 + *r)£ + Y = Y *M{0, 1}) = 1. (20) 

k>2 k -1 

Substituting £ = exp(p) we obtain the discrete fixed point equation (6). For activi¬ 
ties in the fluid domain, there is therefore a unique solution that stays away from the 
radius of convergence, so that the expected bulb lifetime p = (1+Zi)£+ ^ z k^ k 
is finite. Renewal theory then tells us that £ l Sl(z) converges to 1///: the proba¬ 
bility that a light bulb has to be renewed at time L given that there was a renewal 
at time 0 converges to the inverse of the average bulb lifetime. In particular, 
j- log El (z) —>• — log£ and it follows that the pressure is indeed p = — log£. 

For the continuous system we multiply the partition function with exp(— XL) and 
distribute the additional factor exp(— XL) over the rods and empty spaces between 
them. This attributes the weight ^2 k Zk l(r > £fc)exp(—A r) to an interval going 
from the starting point of a rod to the next one (instead of treating emtpy space 
as monomers, we attach an empty interval to the preceding rod). These weights 
define a probability distribution if 

f oo y 

/ Y z fe 1 ( r ^ 4) exp(—Ar)dr = ^Y Zk ex P(~^4) = 1, 

Jo k A k 

and we recognize the fixed point equation (|5|). 
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4. Computation of the pressure and packing fraction 


Here we prove Thcorems l2.2l and l2.3l for the continuous system and Theorem l2.10l 
for the discrete system. 

Proof of Theorem, ! 2. 21 For A > 0, let F( A) := / 0 °° exp(—AL)2i(z)dL be the Laplace 
transform of the grand-canonical partition function with respect to the system 
length L. We note that F( A) < oo for A > p(z) and F( A) = oo for A < p{z). Let 
0(z) be the solution to ^2 k Zk exp(£fc0) = — 9 , if the solution exists, and 6(z) = 9* 
otherwise. 

Ordering rods from left to right and singling out the right end point x of the left¬ 
most rod, we see that the partition function satisfies the renewal equation IFel711 
Chapter XI] 

rL 

5l(z) = 1 + V'zfe / S L _ x (z)da;. 

k •' e * 

It follows that for all A > 0, 

F(A) = i(l + F(A)^z,e- v '“). 

k 

If F( A) < oo, then ^2 k z^e~ Xik must be finite as well, thus A > — 0* and 

(a - y^ Zfc exp(-4A))F(A) = 1. (21) 

k 

Since F( A) > 0, we deduce J2k Zk ex p(—A^*,) < A which implies that —A < 9{z). 
Conversely, if —A < 9(z), then A > ]T] fc Zk exp(—A£fe) and F( A) is the unique 
solution to Eq. m, so in particular F( A) < oo. 

We have shown that F( A) < oo for —A < 6{z) and F( A) = oo for —A > 9(z). It 
follows that p(z) = —0(z). □ 


Proof of Theorem 12.10! For the discrete system, the grand-canonical partition func¬ 
tion satisfies the discrete renewal equation [IVZ061 

L—l 

El(z) = (1 + zi)Hl-i(z) + Y z fc S L _ fc (z) + z L . (22) 

k—2 

Eq. (ED for the Laplace transform in the continuous setting is replaced by a relation 
between generating functions 

OO OO OO 

1 + Y = ((1 + Zi)£ + Y ~fc£ fc ) (i + Y 

L—l k—2 L =1 


which yields 


1 + Y = 


L—l 


i-(i+^-Er = 2 z ^‘ 


The radius of convergence R on the left-hand side is given by the solution R < 
exp(0*) to (1 + zi)I? + 2 z kR k = 1 if the solution exists, or equal to exp(0*) if 

it does not. Since p = — log I?, we are done. □ 


Now we come to the limit law for the packing fraction of the continuous systems. 
In principle, we could deduce it from Theorem 12.31 using //-dependent activities 
Zfc(/i) = Zfcexp(£fc/x) and general theorems relating limit laws to differentiability 
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and strict convexity of the pressure as a function of [i. We prefer to follow a more 
direct approach, which has the advantage of shedding additional light on the value 
of the pressure. In particular, Lemma l4.2l below expresses the pressure in terms of a 
variational problem whose minimizer corresponds to the limiting packing fraction. 

We start by computing the multicanonical partition function. For N = (N k ) k ^ £ 
X* and L > 0, let 


Z l (N u N 2 ,...) 


1 

n k N k '. 


L—ii 

dzn 



dxfcAfc • • • l(rods do not overlap) 


(23) 


so that Sl(z) = 1 + X^ivez* zN % l{N). The next lemma is a variant of the 
representation of the canonical partition function for non-overlapping rods of fixed 
length |Ton36l . 


Lemma 4.1. For every N £ I* and every L > 0 such that N k £ k < L , we have 


Z L (N 1: N 2 ,...) 


{L-T. k NA)^ Nk 


(24) 


Proof. The factorials in Eq. (E51) can be dropped if we integrate only over sectors 
where rods of a given type are labelled from left to right, i.e., x k \ < • • • < x k N k for 
all k. Write M := Y) k Afc- The partition function Zl(N) is a sum of integrals of 
the form 

) rX2—£k(2) 

d xm / & xm -\ ■ ■ ■ / dxi. 

H-Mfc(M) •> i)H-Mfc(M-i) J l) 

Here Xj represent the end points of rods, and the sum is over color assignments 

A;(l),..., k(M) compatible with Ni,N 2 , _ There are (^ ^ ) such assignments, 

and each integral equals M! -1 (L — Nk^k) M ■ □ 



Next, we note that Eq. (1241) can be reinterpreted in terms of Poisson random 
variables. Given S > 0 and 9 £ (—c»,#*), let N k (uj), k £ N, be independent 
Poisson random variables N k (uj ) ~ Poiss((L — S)z k exp (9£ k )) defined on some com¬ 
mon probability space (H, .F,Pfii,s). A straightforward computation shows that if 
N\, N 2 , ... are integers (in No) such that ]T) fc N k £ k = S , then 


z n Z l (N 1 ,N 2 ,...) 

= exp((L - S) J2 Zke eik - 05)p e ,s(vfc € N : N k (u) = N k ). (25) 

k 

Note that the right-hand side depends on 9 , but the left-hand side does not. Eq. (l25l) 
suggests that in the grand-canonical ensemble, the packing fraction satisfies a large 
deviations principle with convex rate function 1(a) + p(z) where 


1(a) 


sup 


9a (1 


6 < 6 * 


o)g(0) 


( 26 ) 


Remember g(9) and 9* from Eqs. (flUl) and (fill) . Let u* := lime /'B* g'(0)- Set 


* 


U* 

1+U* ’ 


1 , 


a 


if u* < oo, 
if u* = oo. 
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Lemma 4.2. We have 

min 1(a) = —p(z). 

<t€[0,1] 

Moreover 

(a) In the fluid regime 1(a) has a unique minimizer a(z) £ (0, a*) given by 
a(z) = Efc4^fcexp(-4p(z))]/[l + X)jfc4-2fcexp(-4p(z))]. 

(b) In the transition regime the minimizers of 1(a) consist precisely of the 
elements of [a *, I]. 

(c) In the close-packing regime 1(a) has the unique minimizer cr = 1. 


Proof. Let <p(u) := sup egR [0u — g(9 )] be the Legendre transform of g(9), where 
we agree g(9) = oo for 0 > 9*. In (0,u*), <p(u) is strictly convex and smooth; 
its derivative 9 = <p'(u) solves g'(9) = u, and as u \ 0, we have ip(u) —> 0 and 
9 = ip'(u) —> —oo. On [it*,oo), ip is affine with slope 9*. We note 


inf I (a) 
<re(o,i) 


inf (1 - a)(p(- - 

<re(o,i) v 1 — a 


) = inf 


ip(u) 


>0 1 + U 


For u £ (0, u*), we have 

d ip(u) ip'(u)(l + u) — tp(u) 
du 1 + u (1 + u) 2 1 ^ 7 ) 

-^-(ip'(u)( 1 + u) - ip(u)) = <p"(u)( 1 + u) > 0. 
du v 7 

Hence the derivative of p(u)/(1 + u) vanishes at u if and only if tp'(u)(l + u) = ip(u). 
Write 9 = ip'(u) and remember g'(u) = 9. Then the equation for a vanishing 
derivative becomes 9u + 9 = ip(u) = 9u — g(9), i.e., 9 = —g(9) and we recognize our 
old fixed point equation. 

In the fluid regime the fixed point equation has the unique solution 9(z) = 
—p(z) < 9*, and u = g'(9(z)) = Ylk 4~fc ex p(~ ^kp(z)) £ (0, u*) is a strict local 
minimizer of p(u)/( 1 + u) in (0,u*). Correspondingly a = u/(l + u) £ (0,<r*) is a 
strict local minimizer of /(•) in (0,cr*). But /(•), as the supremum of a family of 
affine functions, is convex, so a is the unique minimizer of /(•) in all of (0,1), and 


inf I(a') = 1(a) 

cr'e [0,1] 


1 + u 


ip' (u) = 9 = —p. 


This proves part (a) of the lemma. 

In the transition regime a variant of the previous argument shows that 1(a) is 
decreasing on (0,a*) and I(a*) = 9* = —p(z). If a* = 1, we are done. If a* < 1, 
then u* < oo and for u > u* 

<p{u) ip(u*)+9*(u-u*) 9*u* — g(9*) + 9*(u — u*) ^ 

-=-=- =9 = —p(z). 

1+u 1+w 1f« 

It follows that I is constant on [a*, 1]. This proves (b). 

In the close-packing regime the fixed point equation has no solution and / (a) is 

decreasing on (0,cr*). On (cr*, 1) or (u*,oo) we have 

P(n) = -g(9*) + 9*u =&¥+ -g(9*) - 9* 

1 + u 1 + u 1 + u 

which is decreasing in u. Thus 1(a) is decreasing in all of (0,1) and 


1(1) = inf 1(a) = 9* = -p(z). 

o-G(0,1) 
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This proves part (c). 


□ 


Lemma 4.3. Let [< 71 , 02 ] C [0,1]. Then 


lim sup — log 

L—too L/ 


E 

N=(N 1 ,N 2 ,...) 


z N Z L (N)\(^Y. N ^k 

k 


G 



< — min /(ct). 
<Te[CTi,cr 2 ] 


Proof. We first estimate contributions from small intervals [ 0,0 + e] C [ 01 , 02 ]- 
Abbreviate S = N k £ k and remember the measures P e,s defined above Eq. ([25]) . 
Observe 


z N Z L (N) < exp (-9S+(L - La)g(9))v e ,La(yk : N k (co) = TV*) (28) 

for all 9. We choose 9 as the maximizer 9(a) of 9a — (1 — a)g(9). If 9 > 0, we have 
— 9S < —9aL for all S G [ 0 L, (0 + e)L\ and 


Y j z n Z l (N)1 (S/L€ [ 0,0 + e]) 

N 

< expf^L/(a))p^ )La f^y^7V fc (q;)4 g + 

k 

hence 

lim sup — log E z n Z l (N)1(S/L& [ 0,0 + e]) <-I(a). 

L—*oo / N 

If 9 < 0, we have instead — 9S < —9L(a + e) and 

lim sup y log^) z N Z L (N)l(S/L G [ 0,0 + e]) < — 1(a) — 9e. 

L->oo L N 

If 01 > 0, then 9(a) > 9(a\) =: 9\ stays bounded away from — 00 . We split [ 01 , 02 ] 
into m slices of width e = (02 — 01 )/to and obtain the upper bound 

— min 1 (a) + e max(— 61 , 0 ). 
ere [cri ,cr 2 ] 


Letting e —► 0 then yields the desired bound. 

If 01 =0 we have to proceed more carefully. Fix a 0 such that 9q = 9(ao) < 
min(0*, 0). For S G [ 0 , 0 + e] C [0, 00 ] we apply the bound (l28l) to 9 = 6(a + e) < 
6> 0 < 0 and find 

limsup j log^ z N Zl(N)1(S/L G [ 0 , 0 + e]) < — I (a + e) + eg( 6 ). 

L—tco L N 

Note that g(9) < g(9o) is bounded uniformly in 0 G [O, 0 o]. Splitting [0, 00 ] into 
slices of width e and letting first L —> 00 and then e —> 0, we obtain an upper 
bound for the sum over packing fractions in [0, 00 ] of the desired form. This works 
for every sufficiently small 00 . Combined with the bounds for intervals [ 01 , 02 ] with 
01 > 0, this proves the lemma. □ 


Now we can prove Theorem l2.3l 


Proof of Theorem \2.3l By Lemma 22] in both situations (a) and (b) of the theorem, 
1(a) has a unique minimizer 00 , given by a(z) or 0 = 1. Let e > 0. Then by 
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Lemma ll.dl the grand-canonical probability of seeing a packing fraction in [0, ao — s] 
or [(To + e, 1 ] is bounded by 


expf-L min (/(er) + p(z)) + o(L)\ 
V Iff —ffnl>£ ' 


\a—ao\>e 

which goes to zero exponentially fast as L —> oo. 


□ 


5. Computation of the pressure-activity expansion 

In this section we show that the fixed point equation for the pressure determines 
the coefficients of the (formal) activity expansion uniquely. The convergence of the 
expansion is dealt with in Section [ 6 ] Because of alternating signs, it is convenient 
to work with F(z) = —p(—z ). 

5.1. Continuous system. The functional equation © for the pressure translates 
into 

F{z) = Zk ex P(4 F{z)). (29) 

k 

Let F(z) = J2 n a(n)z n be a formal power series. The right-hand side of Eq. (l29l) 
can be read as a formal power series whose coefficients are polynomials of a(n); we 
say that Eq. (l29[) holds in the sense of formal power series if the coefficients are 
equal, i.e., 

VnGl: [ z n ]F(z ) = [; z ri ] ^ z k exp (4 F(z)). 

k 

The notation u [z n ]F(z)” stands for “coefficient of the monomial z^ k in the 
formal power series F(z)”. 

Proposition 5.1. The formal power series F{z) = Y^n a ( n ) zTl satisfies Eq. (l29l) 
if and only if a{n) = ^y(X)fc 4 n fc)^ fc nk ~ 1 for all n £ I*, and a(0) = 0. 

We give two independent proofs, a direct proof based on the fixed point equation 
and Lagrange-Good inversion |Goo60| . and a combinatorial proof based on a con¬ 
nection between Eq. (1291) and recurrence relations between weighted trees. The first 
proof is more robust and is easily adapted to the discrete model. The second proof 
is of interest as it provides yet another example for the deep connections between 
combinatorics and cluster expansions [Far 101 . 

Proof of Provosition \5. 1\ via Lagrange-Good inversion. The identity (l29l) is equiv¬ 
alent to a system of equations for the expansion coefficients a(n). The system is 
nonlinear but has a triangular form, i.e., a(n) is expressed in terms of coefficients 
a(m) corresponding to total degree < Ylk n k- The existence and uniqueness 

of a solution is therefore easily proven by induction over Thus we are left 

with the computation of the coefficients. Clearly a(0) = F(0) = 0. For non-zero 
multi-indices, write Wk(z) := Zkexp(£kF(z)). We have 

z k = w k exp(-4 ^ Wj ), k£ N. 

3 
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We look for the expansion of F(z) = J2j w j( z ) i n powers of the z k s. A formal 
computation suggests 

1 f d z k 


M(E^w)-n(5s 

3 k 

-n(si 


Z * 3 

d w k 


Z k {w) 


C^2 w o( z )) 

5> )de '((£A~ 


n k + l 




The Lagrange-Good inversion formula IGoofiOl IEM94j says that indeed 

OO OO rv 

[*”]<E»*(*» = I™”] de *((sr)« e ,u P pJ}- 


3 = 1 


7 = 1 


The partial derivative is = (5 k t — £ k w k ) exp(—4 ^ t«j) and the determinant 
equals 


det 


(Ske-£kU>k) =1- V 4 Wk- 

\ / fc,^Gsuppn z —' 


fcGsupp n 

Note for matrices A of rank one, det(id + A) = 1 + Xu=i a n ■ It follows that 


[*”KE w i) = ME ^)(i - E 


lUk n ‘ ( ‘)(E, v>j) 


Now 


3 


(E fc "fe4)(E,- »i) _ 


J •" = |tu 


Therefore (|n| = J2 k n k) 




n,(»>jO 


M-i 


(Ej-raj - !)• 

^ Hfc4) 


k 

3 "‘3 1 


;E ®(E^) 


E, ">) 


fe 


Ini - 1 


j?i|-2 


“ r,! (EM 


|n|-l 


□ 


For the combinatorial proof of Proposition 15.11 we identify -F(z) as an expo¬ 
nential generating function of colored labelled weighted trees. First we need some 
definitions. A colored set is a pair (V, c ) consisting of a finite set V and a color 
map c : V —> N. We consider two colored sets as equivalent if for each color k, 
they have the same number n k of elements of a given color; equivalently, if there 
is a color-preserving bijection between them. Given n £ I*, let (V,c) be one rep¬ 
resentative of the equivalence class with color numbers n k - For example, we can 
choose V = {(1,1),..., (l,m), (2,1),..., ( 2 , 712 ), • • ■} with the natural color map¬ 
ping c(fc, n) = k. Let T{n) be the collection of rooted trees on V. We think of 
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such trees as colored rooted trees labelled within colors. To each such tree we assign 
a weight given by 

^(t) • || ^c(father in e) ■ 

e£E(~/) 

Put differently, the weight of an edge originating in a father with color k is £ k , 
and the weight of the tree is the product of the edge weights. The corresponding 
exponential generating function is 


G(z) = £ 

n.£Z* 


Z 

n\ 


Y w< ^)- 

7€T(n) 


Proof of Proposition \5.l\ via combinatorics. The combinatorial proof is in two steps: 
1. show F(z) = G(z), 2. compute the sum of weights ^ 7£7 -( n ) w( 7 ). 

First we note that G(z) satisfies the same functional equation as F. This can 
be seen by summing over the color of the root, and then over the subtrees attached 
to the root, see |Goo65l for similar relations for colored trees. The functional 
equation determines the power series coefficients [z n ]G(z) uniquely, as can be seen 
for example by an induction over J2k nk - As a consequence, F(z) = G(z). 

For the computation of the weights, let di be the degree of a vertex in a given 
graph 7 . Each vertex is the father of di — 1 children, except the root which has di 
children. Thus we can compute the sum of weights by summing first over all possible 
degree distributions and then over the choice of the root. We use the formula for 
the number of labelled (unrooted) trees with a given degree distribution [Moo701 , 
write m = \n\ = n k , and obtain 


E w W 

7 GT(n) 


E 

(di)i<i<m• Ei di=2m—2 


m — 2 

d\ 1 ,..., dm 1 


m m 

(n ^c(i) *) (e ^«) 


i=1 


2=1 


m 

(£■ 

2=1 


m —1 


(E n ^) 5: * n * 1 


Here we have implicitly chosen our trees as trees on V = {1,... ,m}; recall that 
c : V —> N is a fixed map such that #{* £ V \ c{i) = k} = n*,. □ 


5.2. Discrete system. For the discrete system, the fixed point equation after sign 
flip becomes 

OO 

exp (F(z)) = 1 + Zk exp (kF(z)). (30) 

k =1 


Proposition 5.2. 

if and only if 


The formal power series F(z) = a(n)z Tl satisfies Eq. (l30l) 

1 v (E k knk ~l)\ 


a(n) = 


U k n k \ (Efc^fc - Efc«fe) ! ' 


for all n £ I*, and a( 0 ) = 0 . 


We provide a proof based on Lagrange-Good inversion and leave open the combi¬ 
natorial interpretation. Note that G(z) = exp(F(z)) can be interpreted as an ordi¬ 
nary generating function for non-labellcd, colored planar trees. For example, when 
all activities except for a given rod length k vanish, Eq. becomes G = 1 + z k G k 
which is the generating function for k -ary trees (every vertex has exactly k chil¬ 
dren, only non-leaf vertices are counted), and the expansion coefficients of G are 
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generalized Catalan numbers )HP9ll . When k = 2, G = 1 + Z 2 G 2 is the generating 
function of the standard Catalan numbers -Wy ( 2 ”) and is associated with binary 
trees. The interpretation of F = log G, however, is less clear. 

Proof. The identity F = log(l + Efc zu exp (kF)) leads to a system of equations for 
the expansion coefficients a(n) that is triangular, and the existence and uniqueness 
of the solution can be proven by induction just as for Proposition 15.II Thus we are 
left with the computation of the coefficients. 

Let Wk = Zk exp(kF(z)). Then exp (F) = 1 + J2 < kLi w k and Zk = wt/{ 1 + 
Y^jLi w i ) k • We compute 


dz k _ S ki kw k 

dw e ~ (1 + £“ r ~ (1 + r u,j) k+l 

and for every finite non-empty set I C N 

det (—) = _ - _ (l _"l 

The identity a(0) = 0 is obvious. Let n £ 1* be a non-zero multi-index. By 

Lagrange-Good inversion, we have 

\z n ]F{z) 

OO OO rj 

= [w n ] {log(l + 51 ^)(! + 55 «'j) E ‘ £ " p P n Hnk+1) det {(g^)k,jGs u PP ™) } 

0 =1 i=i 3 

= [w n ](i + 55 iog(i +E “i) 

3 3 

- e [^ n - efe ]fc(i+ E w ^ k fcnfc_11 °g ( 1 + E w i)- 

fcGsupp n j j 

Let Hq = 0 and H m = 1 + i + • • • + — be the harmonic numbers. For m £ No set 
f m (u) = ^,(1 + u) m (log(l + u)~ H m ). Then f' m+1 (u) = f m (u) and an induction 
over m leads to 


m , \ 

(1 + u) m log(l + u) = E (J) {Hm - H m _k)u k + 0{u m+1 ). 


Let N = Efc krik and M = J2k nk- It follows that 


[W n ]{l + E Vj) N log(l + E ={ M ) x [“"K 1 + U ) N lo s(l + «) 

} } \ni, n 2 , •. •/ 

M 


N 


(H n — H n - M ) 


ni,n 2 ,.../ \My 
1 AT! 

n! (N-M)'S Hn ~ Hn ~ m ^ 
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Similarly, 


[«»"“' 2 fe ](i+M 1 +E^-) 


M- 1 

, . . . , Tl k 1, . . . 

M — 1 

ni, ...,n k -l,...J \M - 1 
n fc (AT - 1)! 


x ^-^(l + u^-MogCl + u) 
{Hn-i — H n - M) 


Therefore 
[. z n ]F(z) = 


n\ {N-M)\ 
N\ 


(Hn-i — H n -m)- 


n\ (N — M)\ 
N\ 


(Hn ~ H n _ m ) — E - TK\^ N ~ l ~ Hn-m) 


n\ (N — M)\ 


n\(N — M)\ 


(Hn — H n _ i ) — 


(N — 1)! 
n\(N — M)V 


□ 


6. Convergence of the activity expansion 

6.1. Continuous system. Here we prove Theorem 12.41 The fixed point equation 
for p and Proposition 15.11 show that if p(z) has a convergent expansion, then that 
expansion is necessarily given by Eq. m- Thus it remains to investigate the 
convergence of the series. Because of the alternating signs, it is enough to look at 

f(*)= £ f5-(£W) E, ’“" 1 

nGl* k 

for non-negative variables z k > 0. 

The proof of the divergence criterion is based solely on the fixed point equation. 

Proof of Theorem \2.4^ b). Take variables z k > 0. Recall from Section [S] that F(z ) = 
^2 k z k exp (l k F(z)) in the sense of formal power series. If F(z) is finite, the iden¬ 
tity holds as an identity between positive numbers. Since F > 0, it follows that 
inf a (]C fc Zk exp(^fca) — a) < z k exp (£ k F) — F = 0. Taking the converse, we see 
that if inf a (]T fc z k exp (l k a) — a) > 0, then F(z) = oo. □ 

For the proof of the convergence criterion, we use the explicit expression of the 
coefficients and Poisson variables (compare the proof of Theorem 12.31 in Section HJ). 

Proof of the convergence in Theorem Pm (a)- Take z k > 0. For V € V := {J2 k n fc4 \ 
n £ I*}, let 

h v -= E —p^i(e^4 = v). 

n£Z* k 

Note F(z) = J2vev h v/V. 

Assume first that J2 k z k exp(a£ k ) < a for some a > 0. For V > 0, let R k ~ 
Poiss(P 0 feexp(a£fc)), k £ N be independent Poisson random variables defined on 
some common probability space (fi, F, Py). Note that 

h v = e-“^Pv(E 4-Rfc(w) =V), a := a - ^2 z k exp (al k ) > 0. 

k k 
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More generally if V £ V D [n, n + 1) for some n £ N, then 

h v < ^(^ + l) Efcnfc l(E nfc4 = y ) 

nEl* k 

= P n+1 £ k R k (u) = V ) e ~ aV+ ( n+ V E 

k 

< Pn+l (J2 Zk R k(u) = V) e Q_ “( n+1 ). 


/c 


fc 


exp(af fe ) 


Therefore 


CXJ , CXJ ^ 

E E ^<Er" K1, <“ 

n=l V6Vn[n,n+l) n=l 

Note that if £k —t 0, the set V £ Vn[n,n + 1) can be countably infinite. In addition, 

H = e Mj2 Zk ) - 1 < e a - 1 < 00 . 

VGV: V<1 n£l* ' k 


It follows that F(z) < oo. Next we show, still under the assumption Zk exp (a£k) < 
a, that 

F(z) < a. (31) 

To this aim consider the curve y = Ek Zkex p(^fc$), $ > 0. It is the graph of 
a convex increasing function that starts at g( 0 ) > 0, i.e., above the line y = 6 
and is below the line at 6 = a, g{a) < a. It must cross the line at some 0\ < a 
and possibly crosses it again at some $2 > a. Moreover g'{9 1 ) < 1 < g'(6 2 ). We 
already know that F(z) solves g(F) = F. If the fixed point equation has a unique 
solution, then F = d\ and we are done. If the equation has two solutions, we 
invoke a continuity argument: for t £ [0,1], write tz = (tzk)k&n- Then F(tz) solves 
^ z k exp(£fc(?) = 6 . For t close to 1, the latter equation has two solutions 0i(t), 
62 (i) with $ 1 , 2 ( 1 ) = $ 1.2 and 


J2^tz k e e ^ < 1 < Y,£htz k e 6 ^ e \ (32) 

k k 

It follows that as t \ 0, either 0 2 (t) ceases to exist or 9 2 {t) 00 . On the other 

hand the power series expansion for F converges for all t < 1 and has no zero-order 
term, thus 1 1 —>■ F(tz) is continuous on [0,1] and converges to 0 as t \ 0. It follows 
that F(tz ) = 0i(t) for all t £ [0,1]. Therefore F(z) = 0\ < a. This proves Eq. (I3T1) . 

Now consider the case Zk exp (a£k) = a for some a > 0. Let t £ [0,1). Then 
Ek(^ Zk ) exp (a£k) =ta < a, hence 

~71 

F(tz) = 'V' —r(Y^ n k (.k)^‘ knk ~ 1 i^‘ knk < a. 

z ' n\ z —' 

n k 

We let t /* 1, use Abel’s theorem and conclude that F(z) < a < 00 . □ 


6.2. Discrete system. Here we prove Theorem 12.111 Because of the alternating 
signs and Proposition 15.21 it is enough to look at 


\ _ ST' zn 

n! (Ek kn k-Ek nk V- 


for non-negative z k - 
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The proof of the divergence criterion is based on the fixed point equation and is 
completely analogous to the continuous setting. 

Proof of Theorem \2. 1 lV b) . If F(z) < oo, then it satisfies the fixed point equation 
exp(.F(z) = 1+X^fc z k exp(A:F(z)). Setting a = F we find Ek Zk exp(afc) < exp(a) — 
1. It follows that if Ek z k exp (ah) > exp(a) — 1 for all a, then F(z) = oo. □ 

For the sufficiency of the convergence criterion, we have to replace the Poisson 
distribution employed in the continuous setting by something else. 


Proof of the convergence in Theorem \2.1lV a). Assume first that 


for some a > 0. Let 


1 + Ek> l z k exp (ka) 
exp (a) 


exp(— a) < 1 


Po ■= 


1 

1 + Ej>i z j exp(aj) ’ 


Pk ■= 


Zk exp(afc) 

1 + Ej>i exp(aj) 


(k g N). 


Fix S G N. Consider S independent random variables with values in No and 
probability distribution (pk)- We may think of S boxes that are either blank or 
have the color k, and each box chooses its color independently according to the 
probability weights pk- Let (nk)keN be a sequence of non-negative integers with 
M = Ek nk < S. The probability that S — M of the boxes are blank and the 
remaining M boxes have exactly nk boxes of color k, for all k £ N, is given by 



S 

M,m,ri2, ■ ■ ■ 


Po 


n 

k> 1 


G - M, nu n 2 , ...) (1 + Ej z i exp (aj)) s 

S\ z n / exp(a) 

(S - M)\ n\ V1 + Ej Zj exp (aj)J 


It follows that 


E 

«.f=Z* 


Z n (Ek kn k~l)'- 
n\ {Ek kn k-Ek n k) ] - 




exp (—aS) 
S 


and 

S>1 

In order to cover the case that Efc>i z k exp(fca) = exp(a) —1 for some a > 0 but “<” 
fails, one exploits Abel’s theorem and the fixed point equation. The procedure is in 
two steps: 1. Backtrack to the case that Ek z k exp(fca) < exp(a) — 1. Then F(z) is 
finite and satisfies the fixed point equation F = log(l + Efc>i z k exp(fcF)) =: h(F). 
If h(8) = 9 has two solutions, use analyticity to deduce that F(z) is equal to the 
smallest solution and that F(z) < a. 2. When Ek ^feexp(fca) = exp(a) — 1, go to 
(tzk)km and exploit Abel’s theorem for t 1. The details are analogous to the 
continuous setting and we leave them to the reader. □ 
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7. VlRIAL EXPANSION 

Here we prove Thcorem l2.5l for the continuous system. The proof of Theorem l2.12l 
for the discrete system is similar and therefore omitted. A heuristic derivation of 
the relevant expressions for both the continuous and the discrete systems is given 
in Section EH] 


Proof of Theorem, \2.5\ Fix j £ N and z k , k ^ j. The inequality Zk exp(—f k 6*) > 

—0 * extends by continuity to some open neighborhood of Zj. and in this neighbor¬ 
hood the pressure solves the equation F( Z j,p) = 0 with F( Z j,p) = p—^k z k exp(— p£ k ). 
Since 


dF 


— (zj , p) = 1 + z k he > 1 > 0, 


(33) 


the implicit function theorem shows that the partial derivative exists and sat 


isfies 


dp Zj exp(— p(z)£j) 


Pj ( z ) ~j g z . ( z ) 1 + J2 k z k £ k exp(-p(z)4)' 


We multiply with A. sum over j , and get 


r ( z ) = ^Z^pA z ) = 


E.,- Z jtj exp (-p(z)tj) 
zjfj exp{-p{z)£j) 


< 1. 


(34) 


(35) 


If we sum Eq. (1341) right away and use the fixed point equation for p(z), we get 

p( z ) 

k 


EftM = TTE 


k z k t k exp (~p(z)£ k ) 


= P( Z )( 1 -^V-kPktz)), 


which yields the expression of the pressure in terms of the densities. Eq. (1M1) also 
shows 


Zj = Pj(z)e?W* (1 + ^ z k l k e~ p ^ k ) 
k 


Pj( z ) 1,-pQs) 
1 - £ fc 4pfc(*) 


□ 


8. The inverse function theorem does not apply 


Proof of Theorem \2.8\ Suppose by contradiction that p(-) is a bijection from U a 
onto 14 with U a and 14 neighborhoods of the origin in suitable spaces E a , Eb- 
Let e > 0, /ife := ek~ 3 exp (—bkf We choose e > 0 sufficiently small so that 
C := Yh k kpk < 1 and p = {pk)ke n G Vf>. Since p(-) is assumed to be a bijection 
from U a onto 14, there is a unique z £ U a such that p = p(z). Theorem 12.51 shows 
that this activity vector is given by z k = (1 — C)~ 1 p k exp (kp(z)) > 0, and we have 
p(z) > 0. Therefore 


1 

1-C 


5> fc e“ fc 


£ V — P (p( z )+ a ~ b ) k 
1-C ^ k 3 

k 


< oo 


and b > p + a > a. Next set z k = — ek~ 3 exp(— ak), with e > 0 small enough 
so that z £ U a - Then p(z) < 0 and p k {z) is proportional to Zk exp(— kp{z)) = 
£/c -3 exp(— (a + p(z))k). Since p(z) £ Vb we must have b < a + p(z) < a. Thus 
a < b and b < a, contradiction. □ 
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